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ABSTRACT 



Signal Processing for Antisubmarine Warfare is a short 
course in electrical signal processing fundamentals and their 
applications in the field of antisubmarine warfare. It con- 
tains an introduction to Fourier transforms and their proper- 
ties, sampling and quantization, filters and bandwidth 
requirements, random signals and noise, and an introduction to 
four types of processing equipment; the DELTIC, energy detec- 
tors, correlation detectors, and beamformers. Course objectives 
are given in terms of specific questions which a person com- 
pleting the course should be able to answer. The course text 
and illustrative material is contained in the appendix to the 
thesis. The course is designed to be presented in the Fleet 
to the personnel involved with the operation and employment 
of detection equipment to provide them a better understanding 
of the operations accomplished by their equipment and to develop 
in them a better appreciation of the problems and limitations 
associated with signal detection in the antisubmarine warfare 
environment . 
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INTRODUCTION 



Signal Processing for ASW was developed as an introduc- 
tory course intended for the Naval Officer who has had the 
traditional training in ASW hardware operations and ASW ad- 
ministrative requirements. However, it is applicable to any- 
one involved in ASW - officer, enlisted, or civilian - who 
does not have an engineering background. 

There are several fields that are of particular interest 
to the ASW specialist: acoustics, signal processing, environ- 
mental sciences, and probability/statistics. This course 
addresses only the signal processing aspects of ASW. Short 
courses covering tlie uLher three aieas are also under develop- 
ment at the Naval Postgraduate School. Other theses relevant 
to these areas are KING and SWIFT, 1975; DECKER and SOLLEN- 
BERGER, 1975; and SMITH and VICE, 1975. For anyone who takes 
all four of these courses, it is recommended that the signal 
processing be taken last, because the information contained 
in the other three will enhance the student's ability to 
obtain an understanding of the capabilities and limitations 
of signal processing. 
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COURSE OBJECTIVES 



The overall objective of this course is to familiarize 
the student with the principles of signal processing used in 
ASW sensor systems. The student who successfully completes 
this course will be knowledgable in the specific areas which 
follow to the degree indicated. 

A. SIGNAL PROCESSING THEORY 

Understand the theory of signal processing sufficiently 
well to do the following. 

1. Explain the relationship of the Fourier Transform to 
the "time domain" and the "frequency domain". 

2. Given a Fourier Transform pair, identify the opera- 
tions they express. 

3. Given a square wave and a rectangular wave: 

a. sketch the correlation of one with the other; 

b. sketch the convolution of one with the other; 

c. transform the given waves from the time domain 

to the frequency domain, multiply them, and apply the inverse 
transform to return to the time domain; 

d. compare the results obtained in a, b, and c. 

B. SIGNAL PROCESSING FUNDAMENTALS 

Be sufficiently knowledgable of signal processing funda- 
mentals to do the following. 

1. Show the difference between a digital and an analog 
signal . 
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2. Discuss how noise limits quantization in the sampling 
and quantizing of a signal. 

3. Relate Sampling Theorem, Nyquist Rate, and "aliasing”. 

4. Given a plot of signal and noise as power density 
versus frequency, show how filtering increases the ratio of 
signal power to noise power. 

5. Given a plot of filter response versus frequency, 
identify a high pass, a low pass, and a bandpass filter. 

6. Relate filter bandpass, integration time, and 
frequency resolution. 

7. Define the statistical characteristics of a random 
signal and relate them to the signal voltage (or current) 
o-nd p c A* cornponciicS * 

8. List the assumptions necessary to apply statistical 
techniques to the processing of random signals. 

9. Contrast the cross correlation of correlated signals 
with the cross correlation of uncorrelated signals. Show how 
the correlated signals are enhanced. 

10. Explain the three effects of filtering gaussian 
white noise. Discuss the relationship of these effects when 
processing signal and noise. 

C. DELTIC (DElay Time Compressor) 

Be able to describe the DELTIC, including; 

1. how the DELTIC functions; 

2. the benefits of using a DELTIC in the processing of a 
signal ; 
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3. given the specifications of a particular DELTIC, cal 
culate the integration time required and the frequency reso- 
lution of the output. 

D. ENERGY DETECTION 

Describe energy detection by explaining: 

1. how the system determines that a target is present; 

2. how the threshold setting affects the probability of 
dedection, P(d) , and the probability of false alarm, P(fa); 

3. the significance of the direct relationship of a 
P(fa) with a given P(d) as represented by a Receiver Opera- 
ting Characteristics (ROC) curve. 

E. CORRELATION DETECTION 
Be able to: 

1. describe the method used to discriminate between the 
returning echo and noise; 

2. describe the effect of doppler on correlation; 

3. explain how doppler is determined; 

4. describe how pulse length effects doppler resolution 

5. describe how pulse length effects range resolution; 

6. list the advantages of using an FM slide; 

7. state the advantages and disadvantages of using 
Pseudo- Random Noise pulses for search and tracking; 

8. given a tv\?o-hydrop]ione array being used with a corre 
lation detector, sketch and explain how the signal direction 
is obtained, why a direction ambiguity exists, and how the 
ambiguity may be resolved; 
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9. explain the limitations on the type of signals which 
may be processed passively by a correlation detector, and 
why these limitations exist. 

F. BEAM FORMING 

1. Describe the method used by the beam former to obtain 
direction. 

2. Name the type detector used with a passive beam 
former . 

3. Given a hydrophone array in a beam former configura- 
tion, sketch the output of the beam former for a signal in 
the beam and a signal out of the beam. Show by comparing the 
outputs the importance of setting the threshold to determine 
a target present in the beam. 

THE STUDENT 

It is expected that the student will be currently in- 
volved in ASW. Prior courses in probability, college algebra 
and calculus, and basic electricity are helpful but not re- 
quired, since the course is self-contained. 

THE INSTRUCTOR 

There are three basic qualification needed by the instru- 
tor for this course: a background in electrical engineering 
that includes the fundamentals of signal processing; a know- 
ledge of current practices in ASW; the ability to teach. The 
Naval Postgraduate School has faculty available with these 
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qualifications, and therefore should be considered as a logi- 
cal source. NPS graduates may also provide a source of people 
with the requisite qualifications. 

COURSE DEVELOPMENT AND PRESENTATION 

This course was designed to provide a thorough understand- 
ing of the principles involved in the application of spectrum 
analysis to ASW signal processing, and the operation of the 
DELTIC, energy detectors, correlation detectors, and beam, 
formers. These processing methods comprise the core around 
which the U.S.NAVY's ASW sensor systems are built. In pro- 
ceeding toward the above goal, the course follows a logical 
approach of building from the basics of Fourier Transforms 
and Fourier Transform Properties. Following those two sec- 
tions is a section on Sampling and Quantization and a section 
on Filters and Bandwidth Requirements. These latter two 
sections may be interchanged, but both are required for the 
section on Random Signal, Power Spectral Density, and Noise. 
The methods of processing are presented last; in order they 
are, DELTIC, Energy Detection, Correlation Detection, and 
Beam Forming. Where appropriate, specific hardware has been 
mentioned, but it is emphasized that the course is primarily 
devoted to principles and methods that may be applied to var- 
ious systems employed on either air, surface, or subsurface 
platforms. 

This course is designed to be taught in the Fleet. 

Although it requires no outside references, each student 
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should be provided with a complete copy o£ the course materi- 
al for his (her) use and retention. The course objectives 
are very specific. They are designed for use by the instruc- 
tor and student as a guide for measuring achievement: they 
indicate to the student \>rhat he should learn; they should 
serve the instructor as the basis for examination questions. 
The only training aid required by the instructor is a chalk 
board, as students may follow along on the illustrations in 
the course material. It will be helpful if 35 mm slides of 
all the illustrations are available. Overhead projectors, 
although common at shore stations, are not common on board 
ship, whereas a 35 mm slide projector is normally on board 
for other training requirements. References to books which 
covei' the subjects in greater cletaxl, arid vvlixch liave been oT 
assistance to the authors in the preparation of this course 
are included in the bibliography. 

If the course is presented independently, no more than 
two hours per day should be devoted to in-class presentation. 
The student should plan to spend a minimum of one hour of re- 
view and study for each in-class hour. The student must have 
time for the information to become impressed in his memory 
and become familiar to him. If all four courses are present- 
ed concurrently (feasible in a training environment due to 
the availavility of the students time and approximately the 
same time required for each course) , no more than one in- 
class hour per day should be spent on any one course. Regard- 
less of the method used, the instructor should be available 
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at other than schcduled-class time to 
which the student may have questions. 

The course is designed for twenty 
in-class time requirement schedule is 



clarify points on 

in-class hours. The 
as follows. 



SECTION HOURS 

Introduction to Fourier Transforms- -- 4 

Fourier Transform Properties 4 

Sampling and Quantization - 1 

Random Signal, Power Spectral Density, and Noise- - 4 

DELTIC- - - 1 

Energy Detection- ----------------- 2 

Correlation Detection- -------------- 2 

Beam Forming- ------- -- 1 



T 1 - —— ^ T ^ ^ .-X J t-m y-x/— .xx-.-x.'*-."! ^ -v-x 4- 4- X A- Vx •-» -4- /-• /-^ "X r /— \ -V* /"s "t 

lie jiidL.ei±a.-L eUveieQ eeUj-vcij-eiii, ow c-jicxl. -L.il 



a two-credit hour, one quarter course at NPS. The sections 
which required the most time to develop were the two 
sections concerning Fourier transforms and the section 
on Random signal, power spectral density, and noise. 

In these sections the emphasis is on development of theory 
and concept which proved to be the most difficult to 
explain clearly and in a manner to preclude misunderstanding 
The other areas were easier to address as they lend them- 
selves much more readily to a descriptive explanation. 
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CONCLUSIONS 



This course effectively covers the major areas of signal 
processing applicable to ASW. Its development has been di- 
rected toward the understanding of theories and principles 
used in current ASW signal processing methods. These theo- 
ries and principles will remain, although in future proces- 
sing methods the relative importance of their individual 
applications may be altered. For example, the current trend 
in ASW signal processing is toward more utilization of digi- 
tal processing methods. This will increase the importance 
of thoroughly understanding sampling and quantization. 

An area which remains to be addressed is the use of 
digital computers in signal processing methods. This should 
include an introduction to boolean algebra, elementary logic 
and set theory, the use of truth tables and Karnaugh maps, 
binary arithmetic and electronic logic circuits. It is 
recommended that such a course be developed as part of the 
continuing education program at the Naval Postgraduate School. 
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APPENDIX 



I. INTRODUCTION TO FOURIER TRANSFORMS 

An introduction to signal processing requires, for 
understanding of many processing schemes, a working fam- 
iliarity with Fourier transforms. To gain an understanding 
of what Fourier transforms are, what they can do, and 
how to use them, let us first examine a common type of 
signal. The electrical power used in the United States is 
60 Hertz (cycles per second) alternating current. That is, 
the current changes its direction of flow in such a manner 
as to complete a cycle sixty times per second. If we plot 
the amplitude of the current in a given direction as a 
function of time, the result is a sinusoid. 




of time, it is said to be the "time domain" description. 
If we plot the signal as a function of frequency, we get 
the "frequency domain" description. 
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Since the signal consists of a single frequency, there 
is but one frequency represented on the f-domain plot. 

Let us nov; look at a different signal, say an alternat- 
ing signal of ISO Hz. The time domain and frequency domain 
plots are as follows. 







a 

CL 

t 



o ISO 



Note the way in A\rhich the changes in the signal are evidenced, 



both in tlie time domain and the frequency domain. 
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Likewise, two other signals, 300 Hz and 420 Hz. 
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So far we have looked at signals of a single frequency. 



We are going to look at signals made up of many frequency 
components, but first let us review the concept of phase, 
PHASE: The phase of a sinusoidal signal with respect to a 

reference signal is the relationship between corresponding 
parts of their cycles in the time domain. For example, take 
a signal for reference. 



Comparing a signal to this reference we can determine 
the phase of the new signal with respect to the reference. 

If the new signal crosses the t axis at the same time and 
in the same direction as the reference, it is said to be 
"in phase" with the reference, or to have a phase angle of 
0° , If the new signal crosses the t axis on the down-swing 
when the reference signal crosses on the up-swing, as shown 
here, it is said to be "180°out of phase". The reason for 
using "180° will be shown later in this development. For 
the moment, however, just accept the fact that a signal is 
said to go through a phase angle of 360° each complete cycle. 
The phase of a signal can either lead or lag tlie reference. 
That is, the signal may cross the t-axis on the dov^ui-swing 
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before or after the reference does. The first case would 
be leading and the second, lagging with respect to the ref- 
erence. Because of this, a signal is not generally described 
as having a phase angle greater than 180° since a 180° lead- 
ing signal looks just the same as a 180° lagging signal. For 
example, a signal which is leading the reference by 270° 
looks the same as one lagging the reference by 90° , and is 
usually described as the latter. 




Signals with Multiple Frequency Components: Returning to 
the signal with more than one frequency component, let us 
look at the way in which tliis would be represented 
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in the time and frequency domains. Take, for example, the 
four frequencies we looked at earlier: 60, 180, 300, and 
420 Hz. If we had a signal which was composed of equal am- 
plitude components of each of these frequencies, the fre- 
quency domain plot would look like this. 




Now let us look at the time domain plot of this signal, 
starting with just the 60 Hz frequency present and then as 
WG ndd t’.o the judiciously chosen mnplitudo nnd phssG 

components of the other frequencies. First, take the 60 Hz. 




] 8 



Now add one-third the amplitude of the 180 Hz signal but 
shift it 180° out of phase with the 60 Hz signal. 




Note the fact that the 180 Hz component was added 180° out 
of phase shows up on the frequency domain plot as a negative 
amplitude. The reason for this will become apparent shortly. 
For now, however, let us continue with the other frequencies. 
To the signal we now add one-fifth the amplitude of the 300 
Hz component, in phase with the 60 Hz signal. 
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Now add one-seventh the amplitude of the 420 Hz component, 
but 180° out of phase from the 60 Hz signal. 
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I£ v/e continue to add carefully selected frequency components 
to the signal, the time domain signal becomes more and more 
like a square wave. With 20 components, the wave looks like 
this . 




From this, one can see that by careful choice of component 
frequencies, and their respective amplitudes and phases, 
that a square wave can be generated from purely sinusoidal 
components. If an infinite number of components were used, 
the resulting waveform would be a perfect square wave. This 
concept is the basis on which a large part of our signal pro- 
cessing theory rests. There is not time in this course to 
demonstrate it, but it is a basic theorem of Fourier analy- 
sis that any shape v\faveform can be created from combinations 
of pure sinusoidal signals. 

Complex Quantities: In order to simplify the mathematics 

involved later in the course, we must introduce the concept 
of complex quantities. In solving the mathematical equations 
describing waves, one encounters solutions which contain terms 
multiplied by /^. As one can readily see, there is no real 
number \\/hich, when multiplied by itself, will give a number 



21 




I 




4 



▼ w 1 



which is negative. Thus, these quantities are known as 
"imaginary". They are, in actuality, artifacts of the math- 
ematics and are not directly related to tangible quantities, 
however, they prove to be of use in describing the behavior 
of wave phenomena. 

If we plot all "real" quantities on one axis, all "imag- 
inary" quantities on another, and time on a third axis, we 
get a coordinate system like so: 




One of the results of 
a sinusoidal signal in 
the real and time axes 
can be generated by a 
vector rotating with a 
uniform angular veloc- 
ity in the complex 
plane (formed by the 
real and imaginary 
axes) . Looking at 
this vector down the 



the solutions described above is that 
the real plane (the plane formed by 







time axis shows: 
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As this vector rotates it is also moving down the time 
axis, so the path swept out by the tip of the vector de- 
scribes a helical path around the time axis. 



RtKV 
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At this point it should be noted that there is another 
measure of angular rotation knov^;n as the "radian”. There 
are 2it radians in 360°. This measure is used because it 
simplifies the mathematics involved in solving wave equations. 
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Returning to the vector generator, we see that the number 
of times -it makes a complete 360° (2tt radians) sweep in one 
second of time is the frequency of the signal. The rotation- 
al velocity of the vector, m, is the number of radians swept 
out by the vector per second. Thus o) = 2iTf since it sweeps 
out 2 it radians each 360° rotation and makes f rotations per 
second. Thus m is known as the "radian frequency" of the 
signal . 

If we look at the projections of the vector on the real 
and imaginary planes as it rotates, we find that the projec- 
tion on the real plane is a cosine wave and the projection 
on the imaginary plane is a sine wave. Also, we see that 
the vector rotates 360° during one cycle of the cosine wave 
on the real plane. This is where the use of 360° of phase 
angle arose. 







The vector has components in the real and imaginary 
planes. A mathematician, by the name of Euler, has sho\^^n 
the relationship: 

± i 0 . . 

e •' = cos 6 ± 3 sin 6, 
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where j is taken to mean "in the imaginary plane". 

By vector addition, the vector generating our signal is 
the vector sum of the real component, cos o)t, and the imag- 
inary component, j sin tot. 

Vector = cos tot + j sin tot 

We see that this is of the form of Euler's identity if 6 = cot. 
Thus, the mathematical description of the vector is e“^ 

This vector is known as a "phasor". 

We have seen that cos tot is the projection of the vector 
tip path on the real plane. Solving Euler's identity for 
cos cot, we get: 

Cos cot = Real Part of [e“^ “^3 
= Re[e-^ 

. ±i(Ot . . . . 

line convention or using e ' tor cos cot is used in 
EE because it suggests several properties of the cosine 
not apparent in the written notation "cosine cot".] 

What does this mean? The ± sign indicates that there are 
two solutions. This corresponds to two phasors rotating in 
the +C0 and -co directions, respectively. 
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These phasors unwind in opposite directions along 
t-axis . 




Looking at the projection of these phasors on the rea] plane, 
we see coswt. 
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Solving Euler's identity for Cos cot, we obtain 



Cos cot 



gj ^ g-j cot 

2 



We can see this , as the two contributions of the phasors 
always add in the same direction on the real plane, there- 
fore, the + sign in the formula. Looking at the projection 
of these phasors on the imaginary plane, we see sin cot. 







Solving Euler's identity for Sin cot, we obtain 

^jcot ^-jcot 
e*^ — e 



Sin cot = 



2j 



Because the two projections are always in opposite direc- 
tions on the imaginary plane, we get the - sign in the 
formula . 

We have shown how a single frequency plots as a single 
line in the frequency domain. This relates to a line for 
each phasor. But we have seen that the mathematical solu- 
tions establish two phasors for each frequency, one rotat- 
ing at and one rotating at -m . Thus, in actuality, we 
have been showing only one side of the spectrum, the posi- 
tive frequency side. The full spectrum of cos cot is then: 



28 











-w < 





The lines are both on the same side of the ^ axis because of 
the + sign in the formula 

j mt ^ -j wt 

Cos Wt = 

2 ♦ 

The full spectrum of sin wt is. 



-UJ 








< 





Note that the lines are on opposite sides of the tj axis due 
to the — sign in the formula 



Sin cot = 



j cot 



- j cot 



2i 



Note also that the lines correspond to projections on the 
real plane for cos mt and the imaginary plane for sin cot. 
We have previously shown that arbitrary waveforms can 
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be made from combinations of sines and cosines of the proper 
amplitudes, phases, and frequencies. How do we determine 
which frequencies and what amplitudes and phases to use? In 
effect, we correlate the waveform with various frequencies. 

We compare the waveform to a given frequency and see how well 
they match. One measure of this is the point by point pro- 
duct of the two averaged over a period. This measure, known 
as (n^ frequency coefficient) is the following: 



n 




x(t) dt 



where x(t) is the waveform of interest, e is -the cosine 

boipg c 8 3 .nd the int'?grHl serves to sum the pro- 

duct over the period, , of the waveform. This determines how 

o 

closely the waveform matches a given frequency. In order to 
find the components of a given waveform, we must compare the 
waveform to all frequencies. The "spectrum” of the waveform 
is then the summation of C^^'s multiplied by the frequency 
component over all frequencies. 

00 

X(f) e^"nt 

n= -oo 

Using this technique, to analyze a square wave train, we get 
the spectrum as shown on the following page. 
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We find that all the equal zero except for a few. These 

turn out to be harmonics of the lowest component. This is 
because the signal is compared to sines and cosines. If the 
signal is compared to other forms of signals, the C^'s will 
not necessarily be harmonics. This technique of analysis of 

-v-i-rr* <-• TT«-*r» V\-\t 

• «CtVV.^i-N^XJUO »VCtO CAHVA . t 

as Fourier Analysis. The fact that some of the C^^'s are 
negative indicates that, the frequency corresponding to that 
is to be 180 ° out of phase with the positive sign fre- 
quencies . 

Let us look at some square wave pulse trains and analyze 
their spectra. 






\ 

If we solve for the C^'s and plot we get the curve as shown 
on the following page. 
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I£ the period of the pulse train is increased but the length 
of each individual pulse is kept constant, the following 
curves are obtained. 
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Note that the point where the envelope of the C^'s crosses 
zero remains the same, but the number of C 's increases. 

That is, there are more frequency components present in this 
signal than there were in the. first. Let us now increase the 
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period even more, still keeping the pulse length constant. 
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We see that the zero crossings of the envelope remain con- 
stant, but the number of C^^'s increases again. If we \\rere 
to increase the period to the point where there was only one 
square pulse left, that is, T we would find that the 

number of C^'s would be infinite, but the zero crossings of 
the envelope would remain the same. 
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From these examples we can see that the zero crossings 
interval ’of the envelope is related to the inverse of the 
pulse width, and the interval between C^^*s is related to 
the inverse of the period of the signal. This, and other 
properties of the Fourier relations, will be discussed in 
greater detail later in the course. 

The form of the envelope of the C^'s for the square wave 
appears so frequently that it has been given a name of its 
own. The form is (sin x)/x, which is called "sine” (pro- 
nounced "sin-cee”) . 

Fourier Transforms: We have shown that as the period 

of a signal becomes infinite, its number of spectral compo- 
nents also becomes infinite. In this case, the spectrum of 
a signal can be characterized by the formula of the envelope 
of the ' s . Techniques have been developed for finding the 
formula of the envelope without solving for the C^'s and one 
method of doing this is known as the Fourier Transform, de- 
fined thusly: 

CO 

X(f) = ^ x(t) e’j dt 

- 00 

Note that this is quite similar to the formula used for find- 
ing the C^'s. Instead of giving us the individual hoiv- 

ever, this gives us the formula for the envelope of the C's 
in the f-domain. 

If we compute the transform of a single square wave using 
this formula, we get the following result for the form of the 
solution : 
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X(£) 



sin X 

X 



-*■ sine X 



This is the same result we get if we compute the C^'s for 
a square wave of T and solve for the formula of the 

envelope . 

Since we can find the spectrum, or f-domain characteris 
tics of a signal, given its time-domain formula, it is only 
reasonable that we should be able to do the reverse opera- 
tion. This turns out to be the case, and the operation is 
accomplished by using the "inverse transform". 



By substitution in the f-domain formula of a signal and 
then integrating, one can get the time domain characteris- 
tics of the signal. The two transform formulas, forward and 
inverse, are called a "transform pair". One finds the f-do 
main characteristics given the t-domain formula, and the 
other finds the t-domain characteristics given the f-domain 
formula. Later we will see that these are powerful tools 
used in signal processing. 

For the moment, let us return to the cosine wave which 
we looked at earlier and see what happens if we introduce 
modulation. That is, superimpose it on a carrier frequency 
cosine wave similar to AM radio transmissions, 

A modulated cosine wave looks as pictured on the follow- 
ing page in the time domain: 



x(t) 
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where the envelope varies as the cosine wave we are looking 
at and the individual peaks are at the carrier frequency. 
What does this signal look like in the £-domain? We know 
what the f -domain plot of our cosine wave looks like. 




The carrier is just a cosine wave of a different frequency 
so its plot looks like this. 
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Note that both f -domain plots were centered about f = 0. 
Now, if we put the two together, i.e., modulate the carrier 
v/ith our cosine wave, we get the following: 




We find that the resulting signal is of the form Cosine 
(f^ ± f) . The f-domain plot is the cosine wave plot cen- 
tered about the carrier frequency values instead of about 
f = 0. Thus, modulation in the time domain is equivalent 
to a translation of the spectrum in the f-domain from cen- 
tered about the f= 0 axis to centered about the carrier fre- 
quency. Let us look at a modulated square wave pulse train. 




We knov^^ that the f-domain plot of a square wave is a 
sine function. We have just seen that modulation of a sig- 
nal translates the f-domain plot of the signal and centers 
it about the carrier frequency. Thus the f-domain plot of 
a modulated square pulse train is as shown on the following 
page. 
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where the envelope takes the shape of the square wave spec- 
trum but is centered about the carrier frequency instead of 
about f = 0. 
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II. FOURIER TRANSFORM PROPERTIES 



• The properties of Fourier Transforms are stated most suc- 
cinctly in the accompanying Fourier Transform Theorem.s. Know- 
ledge of the theorems is important in the interpretation of 
spectra as the theorems express the relationship between 
time-domain and frequency- domain operations. 

In the introduction to Fourier Transforms it was shown that 
the frequency spectrum, X(f), of a signal, x(t), was the Four- 
ier Transform of the signal. That is, X(f) =j^[x(t)] and in- 
versely, x(t )=3 ^[X(f)]. These relations are called "transform 
pairs" and are more conveniently denoted by x(t) ‘‘->X(f) . The 
Fourier Transform Theorems describe the properties of various 
Fourier Transform pairs. 

A. LINEARITY (OR SUPERPOSITION) THEOREM 



If a signal, Xj^(t), is multiplied by a constant, a, 
(as would happen if an amplifier with gain=a were 
placed in a circuit) and a signal, X2(t), is multi- 
plied a constant, b, then their respective line spec- 
tra will also be multiplied by the constant a or b. 




aXj^(t) + bx-,(t) <-fraX^(f) + bX2(f) 
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Linear combinations in the time-domain become linear combina 



tions in the frequency- domain . 



B. TIME DELAY THEOREM 



Given the transform pair x(t) ^->X(f) if 
the signal, x(t) , is delayed by t^ sec- 
onds, as happens when the signal passes 
through a capacitive or inductive element, 
to become a new signal, x(t-t^), then the 
spectrum is modified by a frequency de- 
pendent phase shift to become X(f)e 



The translation of a signal in time changes the phase of the 
spectrum but does not alter the complex amplitude. 

To illustrate the Linearity and Time Delay Theorems, con 
sider a signal V'/hich is a linear combination of delayed sig- 
nals: 



x(t-to)f^X(f)e'^“^o 







y(t) = 2x(t-to) - 4x(t-3to) 



Its spectrum is easily written as: 



Y(f) = 2X(f)e’^“'^o - 4X(f)e’^“'^^o 
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G. 



SCALE CHANGE THEOREM 



It has been shown that a translation of the time origin 
may be accomplished using the Time Delay Theorem. The time 
axis may also be expanded, compressed, or reversed by an opera- 
tion known as "scale change". If a signal, x(t), becomes a 
new signal, x(at) , then: 

(1) x(at) is compressed if a is positive, greater than one 



x(t) = 



a > 1 X (at)= 






tzo 






(2) x(at) is expanded if a is positive, less than one-^ 
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(3) x(at) is reversed in time if a is negative. 
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In the Fourier Transform pair, 

x(at) < > 1 ^1 X(|), 

|a| means the absolute value of a, that is, the numerical value 
only is used, the sign is dropped. 

The scale change theorem expresses the property of recip- 
rocal spreading. If the signal is compressed in time by the 



factor a, its spectrum 
the signal is reversed 
X(-^) X(-£) , which i 
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is expanded in frequency by 1/a. If 
in time (a = -1) its spectrum is 
ndicates a reversal in phase. 





The expansion or compression in the time axis of a signal 
occurs most often in the playback of recorded signals. 

D. FREQUENCY TRANSLATION (MODULATION) THEOREM 

If a signal, x(t) , is multiplied by 

its spectrum, X(f) , will be translated in 

frequency by ±f^. Mathematically stated: 

xCt)e^“ct^X(f-f^), 
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It is more common to multiply a time signal, x(t) , by the real 
part of i.e., by cos cj^-t. From previous work with 

Euler's theorem, it is recalled that 

x(t) cos o)ct = l/2x(t) , 

Taking this result and applying it to the Frequency Translation 
Theorem yields the Frequency Modulation Theorem: 



x(t) cos (i)^t 1/ 2X (£- f + l/2X(£+£j,) 
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Or, if a signal x(t) is multiplied by cos w^^t , its spectrum 
X(£) is shifted up and down in frequency by an amount f(~. This 
is the same result notxccd earlxcr when an Hi' pulsw uxaxix wd^ 
formed by multiplying a rectangular pulse train by cos w^t. 

As noted then, and confirmed by the modulation theorem, time 
domain multiplication becomes translation in the frequency 
domain. As an example,. a rectangular pulse of amplitude A, 
period t, centered at t = 0 and duration x is x(t) = An(t/x) . 




The spectrum is recognized as a SING function, or 



X(£) = A SING ft. 
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Now, multiply that single pulse by cos m^^t : 



x(t) = An(t/x) cos w^t . 



By the modulation theorem, the spectrum is: 




X(f) = 1/2 A SING + 1/2 A SING 




This is the result that is used extensively in single- sideband 
and double - s ideband HF communications. 

E. TIME DERIVATIVE THEOREM 



T ^ 4-Vw 
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signal x(t), the result is multipli- 
cation in the spectrum by jm: 



A x(t)^ > jmX(£). 
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Therefore, differentiation enhances the high frequency compo- 
nents in the spectrum. 



TIME INTEGRATION THEOREM 



Taking the integral 






f a signal x(t ) 



is equivalent to multiplying by (di- 
viding by jw) in the frequency domain: 
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x(t’) dtW ^ X(f) 




Integration, then, suppresses the high frequency components 
but the low frequency components are unaffected. 

G. Convolution 

Take two signals in the time domain t', for example, x(t') 
and y(t ). Now ''flip” the y(t*) signal to the other side of 
the axis and it becomes y(-t') as illustrated below. 



In order to make a comparison between the two signals, multi- 
ply them together at every point in time t' from -<» to +«>. 
Mathematically, this can be expressed as 



z 




x(t') y(-t')dt’ 
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In the above case, anywhere in the t' domain where x(t') 
exists, y(-t') = 0, and where yC“t') exists, x(t') = 0. There- 
fore, their product is zero. 

Since the signals are separated by time t, if a time delay 
of t is introduced making y(-t') become y(t-t'), the signals 
will overlap. If t, the time delay, is varied from to +“, 
the y(t-t’) signal will "slip" across the signal x(t'). The 
result is expressed as a function of time, t; 

p+oo 

z(t) = j x(t’) yCt-t')dt' , 

The above integral is known as the "convolution integral" and 
it expresses convolution mathematically, i.e., the area of 
the product of the two signals. A more convenient notation 
for convolution is: 



z(t) = x(t) * y(t) , 

A graphical interpretation of convolution, as shown on the 
following page, may help to see what happens to the signals. 
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The shaded area on the 
time plot of both signals 
represents the area of their 
"common product" for a partic- 
ular value of time delay, t. 
The sliaded area on the z(t) 
plot represents the summation 
of the area of their "common 
product" as the "flipped" 
signal, y(-t), is "slipped" 
across x ( t) . 



"si ipped" 




Z(t)^ 



Ae% 
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The convolution, i.e., the "flip” and "slip”, of two rec- 
tangular pulses produces a trapezoidal pulse with a base of 
length T, + T 2 » Convolution of two identical rectangular 
pulses produces a triangular pulse with base = 2 t and height = 



A r— ^ 







<-T-* 




Now that convolution in the time domain has been explained 
how does it relate to the frequency domain? There are two con 
volution theorems which relate the two domains. 

Xi.cJ it • i^fj 

Convolution in the time domain trans- 
poses to multiplication in the fre- 
quency domain. 

x(t) . y(t)4->X(f) * Y(f) 

Multiplication in the time domain 
transposes to convolution in the fre- 
quency domain. 



For example, look at the now familiar square pulse shown 



on the following page. 





H. Correlation 

In statistics, when it is desired to know how closely 
one distribution resembles another, the two are overlayed and 
the areas they have in common are multiplied. This "common 
product" represents the degree of similarity, or how well one 
distribution correlates to the other. In signal processing it 
is often necessary to compare one signal with the same signal 
which occurs later in time or with another signal. Signals 
vary in time in contrast to distributions which are constant 
in time. Therefore, in order to completely overlay one sig- 
nal on another, displace the signal in time by x units, then 
vary x to "slip" the displaced signal across the one with 
which it is desired to correlate. With each different x used, 
the "common product" will be a different value. In order to 
obtain one number to represent the value of correlation, sum 
the "common product" over all values of delay, t, and aver- 
age it by dividing by time T (or multiplying by . This 
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defines the "correlation function" R(t) and is described math- 
ematically as: 



R(t) 



lim 
T->’0° T 




x(t+x)dt . 



This is known as a "time average" and is denoted by 

R(t) = <(x(t) x(t+T)^^ 

where the ^ (brackets) indicate ^ x(t + x)dt. 

If the signal is periodic in it need only be averaged over 
one period, , instead of all time. The correlation function 
then becomes: 

^^(t) = It ^(t) x(t+x)dt • 

^ o ^ ^ o 

R(x) is more precisely referred to as the autocorrelation func- 
tion, that is, the correlatioii of one signal with itself. It 
may be desirable, as mentioned earlier, to correlate one sig- 
nal, x(t), with another signal, y(t) , displaced x units in 
time to become y(t+x). This results in the cross correlation 
function, R (x) = <^x(t) y(t + x)% . The autocorrelation func- 
tion is more commonly used. Where "correlation function" is 
used it is assumed to be the "autocorrelation function". 

For a graphical representation, the familiar square pulse 
is used. This, as shown on the following page, represents the 
result of the summation of the "common product" as the delayed 
pulse is "slipped" across the other pulse. This initial part 
of correlation is like the previously described convolution 
except that the signal is not flipped. 



^0 




The next step is to take the time average of the summation of 
the common products. It is seen that, as in convolution, cor- 
relation increases the width of the function by + T 2 , or, 
in the case of a square pulse, from T to 2T. 




To time average, take the 
area = 1/2 bh, or l/2(2T)b^, 
and divide by time, 2T^ 

l/2(2T)b^ _ b’ 

2T 2 • 




What 
;t us see on 



does the correlation 
the following page. 



of 



a 



s inusoidal 



signal 



look like? 






When both signals overlay 
completely, the value is 
large and positive. 



4b 
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A small shift will give 
some negative area but the 
result is still positive. 




A 90° shift makes the 
positive and negative areas 
equal and the resultant" val - 
ue zero. 







A shift of 180° results in 
the entire area being nega- 
tive . 




At 270° the positive and 
negative areas are again 
equal and the resultant 
value is zero. 




At 360°, a complete one cycle 
shift, the signal is back in 
phase and the result is a 
maximum positive value. 
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The correlation function of a sinusoidal signal is another 
sinusoid, and, if the time average is taken, has an ampli- 
tude of 

Up to now there have appeared many similarities between 
convolution and correlation. The graphical interpretation 
has also shown this similarity but much of that is due to the 
choice of symmetrical signals for illustration. Now, look at 
the side by side convolution and correlation of a saw-tooth 
signal . 




Convolution 
"flip" then "slip" 



Correlation 
"slip" only 
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The final results are indicated in the next diagram. 




When "time Averaged” 
becomes : 







1 




'f 



So it is seen that if an asymetrical signal is used, the final 
result obtained by convolution may be quite different from the 
final result obtained by correlation. 



III. SAMPLING AND QUANTIZATION 



In the real world, most of the signals of interest to an 
operator of a detection system are of a continuous, or analog, 
nature. On the other hand, most of the more sophisticated 
processors in use today are digital in operation, and thus can 
not use analog inputs directly. Some means of digitizing the 
analog signals must be used. 

To digitize a signal, one must describe the value of the 
signal at any given time by a digital number. This requires a 
measurement and conversion of the signal value to digital units. 
The question of how accurately we must measure the value of 
the signal in converting it to digital form must be considered. 
Take tlit: following signal, for example: 




This signal has a maximum value of 2.5 and a minimum val- 
ue of 0. Since we must convert the continuous values of the 
signal to a set of discrete values, we must then decide how. 
many discrete values we need in the range from 0 to 3 to de- 
scribe the signal adequately. If we decide to use four levels, 
that is, digitize to the nearest integer value, the signal 



5b' 
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will have the following appearance. 




If we double the number of values and digitize to the near- 
est half integer value, the digitized signal looks like this. 




Note that as the number of values used in digitizing the 
signal increases, the digitized curve more and more closely 
approximates the original signal. The values used in digi- 
tizing the signal are known as quantization levels. In the 
case of a noise-free signal, one can digitize the signal to 
any degree of accuracy by increasing the number of quantiza- 
tion levels. 
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The case of a signal with noise is different, however. 

In this case, there is a random fluctuation of the signal 
value because of noise even when the value represented by 
the signal itself is unchanging. The range of these noise 
fluctuations determines how accurately we can measure the 
value represented by the signal even if we have a perfect 
measuring device. To illustrate this, consider first a sig- 
nal without noise and then one with noise. 



V 


\ 






0 





without HOldE 




We observe that in the noisy signal, the value of the 
signal fluctuates about some mean, or average, value within 
some limits. If the noise is truly random in nature, the 
mean value about which it fluctuates is then the actual 
value of the signal being represented. At any given instant, 
any value of the signal between the upper and lower limits 
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of the noise fluctuation could stand for the actual value 
being represented. The accuracy with which this signal 
can be measured is then limited to the range of fluctuation 
due to noise. For example, suppose the signal of interest 
is a constant 3.0 volts, but because of noise the signal 
fluctuates between 3.5 and 2.5 volts. The most accuracy 
with which we can measure this signal is to the nearest 
integer volt. 

Once we determine how many quantization levels are 
necessary, we mjst then determine how often we need to sam- 
ple the signal in order to reconstruct it accrately. 

Obviously, if the signal never changes value, but remains 
constant, only one sample is needed to reconstruct the sig- 
nal for all times. If, however, the signal is not constant, 
we face a real problem. 

How often do we have to sample a given signal? Well, 
v;e know that one sample is not enough if the signal is not 
constant, and we also know that if the signal is sampled 
continuously (i.e., in an analog fashion) we can reproduce 
it exactly. To sample a signal at extremely short intervals 
is very costly and difficult to accomplish. What we are 
looking for is a compromise, a sampling rate high enough 
that the signal can be reproduced from the samples but not 
any higher than necessary in order to keep cost and equipment 
comp le xi t y down . 

We can see intuitively that, there must be some relationship 
between the rate of change of the signal and the sampling rate 
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necessary to be able to reproduce it. If the signal changes 
very slowly, we need sample it only infrequently. If it 
changes rapidly, however, we must sample it frequently enough 
that it does not have time to make several changes in the inter- 
val between samples or we lose these changes and are thus unable 
to reproduce the signal. If the sample intervals are too long, 
the signal is said to be "undersampled" . This condition, in 
addition to causing us to lose some of the information in the 
signal, can also cause other complications. Perhaps you have 
noticed that in some of the old cowboy movies the stagecoach 
wheels appear to turn backwards at times. This is a case of 
undersampling which illustrates one of the other effects called 
"aliasing”. To explain this, let us look at a sine wave. 
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If we sample this at a frequency less than that necessary for 
reproduction of the signal, we can get a series of samples which 
give the impression of a sine wave of lower frequency. 




59 



This is a problem which must be addressed when analyzing un- 
known signals for frequency components since aliasing can give 
spurious results if the analyzer sampling rate is not high 
enough. This will be discussed further in a minute. In the 
meantime, let us look at the sampling process. 

In the time domain, we can visualize the sampling process 
ideally as the multiplication of the signal by a function known 
as the "ideal sampling wave", which is identically equal to 
zero except at the sample times ®^*^*^* 
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The duration of these spikes, known as "delta functions", is 

infinitesimal, approaching zero, while their amplitude approaches 

infinity, and they are defined such that their area (0 x «>) is 

equal to 1. When the signal of interest is multiplied by the 

ideal sampling wave, we get the product of their areas at each 

point T. Since the area of each delta function is 1, and its 

• • 

width is 0, the product of areas is equal to the value of the 
signal at that time. Thus the products of the two functions 
appear .' as below. 
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Let us look at the frequency domain to see more clearly 
what happens when a signal is sampled. Assume that the signal 
to be sampled is strictly "bandl imited” , that is, it has no 
frequency components outside of some limit. 




2W is the bandwidth of this signal. Although the mathematics 
is rather involved, it can now be shown that the ideal sam- 
pling wave spectrum appears as in the next diagram. 









•s t- 
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We know from the modulation theorem that modulation of a func- 
tion causes its spectrum to be shifted and centered about the 
modulation frequency and otherwise does not change it. Thus, 
the spectrum of the sampled signal appears as follows. 




From this we can see that for strictly bandlimited signals, if 
fg > 2W, there is no overlapping of the spectrum components 
and so the spectrum of the original signal is exactly repro- 
duced. If the sampled signal is then passed through a low-pass 
filter v'hich nas the property of allowing only those frequencies 
below its cutoff to pass, we can exactly recover the original 
signal. 
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It follows from the above that the sampling frequency must be 
> 2W, i.e.-, greater than, or equal to, two times the highest 
frequency component of the signal if we are able to reconstruct 
the signal from the samples. If it is not, the spectrum of 
the sampled signal has overlaps which introduce errors into . 
the reconstruction of the signal. 





These overlaps are the cause of "aliasing” and other incon- 
sistencies in the reconstructed signal. This result is formal- 
ized in the Sampling Theorem, or theorem, of uniform sampling: 

If a signal contains no frequency com- 
ponents for |f| > W, it is completely 
described by instantaneous sample values 
uniformly spaced in time with period 

T < 1/2W. 
s - 

This rate of sampling, f^ = 2W, is known as the Nyquist Rate. 
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It is the absolute minimum sampling rate from which a signal 
can be reconstructed under ideal conditions. 

Until now we have considered only ideal sampling and recon- 
struction. In practice, the ideal is rarely, if ever, realized. 
Practical sampling differs from ideal sampling in three obvious 
respects : 

(1) Practical sampling waves are not composed of delta 
functions but have a finite width;. 

(2) Practical filters are not ideal. 

(3) Real world signals normally sampled are not strictly 
bandl imited . 

In practice these differences are of small enough significance 
that for all normally encountered cases only the last is of 
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IV. FILTERS AND BANDWIDTH REQUIREMENTS 



The ocean is an extremely noisy environment in which 
to search for submarines. Noise from waves, wind, breaking 
surf, shipping, sonic mammals, fish and crustaceans, 
seismic activity, rain, etc., is present to some degree 
at all times. How do we go about sorting out these dif- 
ferent sources and detecting submarines? One way is to 
use filters to eliminate as much of the noise as possible. 

In general, there are three classes of filters relative 
to spectral response: low-pass, high-pass, and band- 

pass filters. Ideal lov/-pass filters allow frequencies 
below their rntoff frequency to pass, and screen out 
all above cutoff. Their spectral response is shown in 
the following diagram. 



High-pass filters allow all frequencies above their 
cutoff to pass and screen out all below. 
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Band-pass filters allow all frequencies between thei 
lov\rer and upper cutoff frequencies to pass and screen 
out all others. 




Filters are characterized by their upper and lower 
cutoff frequencies, or their center frequency and band- 
pass (or bandwidth) , B. The bandwidth is defined as the 
passband width measured in positive frequency only, 
that is, the difference f^ - f^^ = 



B. 
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Of course, no real filters have the sharp cutoff 
characteristics of the ideal filter, so it is customary 
to define the bandwidth of a real filter as the band- 
width between the points v;here the response has dropped 
3 dB from the maximum. These points are known as the 

^ Vi. 4.^ • 




A signal which is passed through a band-pass filter 
appears as in the next diagram. 
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As indicated earlier, filters eliminate as much of 
the noise as possible. How does the filter actually 
improve detection capability? To answer this question, 
we examine the way filters lielp improve detection in a 
large class of detectors, known as energy detectors. 

We discuss the way these operate in greater detail a 
little later, but the general principle of operation 
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is that the detector senses the total power (or energy) 
present in the input, and compares this to a "threshold" 
value. If the input exceeds threshold, the detector 
indicates the presence of a signal, and if the input does 
not exceed threshold, it indicates no signal present. 

In order to see better how filtering helps to detect 
signals in noise, consider the "power spectral density" 
plot of an input. The power spectral density is the power 
in a signal that is carried by each frequency component. 




The total power in a signal then corresponds to the 
integral of the power spectral density over all frequencies. 
The power spectral density thus gives an indication of 
which frequencies carry most of the power in a signal, 
and the relative amount of power at a given frequency 
compared to that at others. The power in a given freq- 
uency band is then proportional to the area under the 
power spectral density curve between the lower and upper 
limits of the band. Comparison of the power present 
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in a signal to the power of the background noise present 
indicates 'the detectability of the signal. Consider 
a signal against a background of "white noise" (white 
noise is defined as having equal amounts of all spectral 
components) to see the gains from filtering. 
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We can see that if the detector examines the entire 
frequency spectrum and senses the total power present 
to compare to a threshold value, the presence or absence 
of the signal does not affect the input significantly. 

On the other hand, if we filter the input and only allow 
a narrow band of frequencies about the frequency of interest 
as the input to the detector, then the presence of a signal 
will substantially change the amount of power in the input. 
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This, in effect, raises the signal- to-noise ratio, 
and thus allows us to detect weaker signals than we could 
detect before filtering. In the case of non-white noise, 
the improvement in signal- to-noise (S/N) ratio can be 
even greater. 




Ambient noise in the 
white noise, and careful 
probability immensely. 



ocean is an example of 
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Energy detectors integrate the energy in the input 
over some time period, some averaging the energy and 
some not averaging. The value obtained in either case 
is compared to the threshold, which is set at the value 
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corresponding to ambient noise when the signal of interest 
is absent. Therefore, one must determine how long to 
integrate the input before comparing to threshold. How 
is this determination accomplished? We saw in the section 
on Fourier transforms that there was a relationship 
between the period of a signal and the separation between 
the frequency components that make up the signal, namely, 
the separation between components is inversely related 
to the period of the signal. 





Thus if we assume that the integration time we choose 
is the period of the input, and that the input exactly 
repeats itself every period ad infinitum , then the minimum 
resolvable frequency difference between components of the 
input is inversely related to the integration time. The 
longer the input is integrated, the more accurately the 
frequency components present are determined. If we 
desire to resolve frequencies 1 Hz apart, we must integrate 
for at least 1 second, for frequencies 0.1 Hz apart, at 
least 10 seconds. 
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Another consideration in detection systera design 
is the fact that most common filters in military appli- 
cations today are composed of reactive elements , such 
as RC/LC/RLC circuits, crystals, or transducer elements, 
which have finite time constants associated with them. 

What does this mean? It means that these systems are 
designed to resonate at a particular frequency or over 
a band of frequencies, and to block or short out others. 
These reactive elements require a finite time to "ring 
up", or establish resonance. Thus, systems do not respond 
immediately to inputs, and this must be considered in 
the system design. The response of these circuits has 
the following appearance. 




The time required for the response to rise to 63 % 
of the maximum is known as the time constant of the system. 
The narrower the resonance, that is, the narrower the 
filter bandpass, the longer the time constant. Thus, 
if we desire to resolve frequencies to 0.1 Hz, we must 
allow enough integration time for the filter to respond 
to each frequency, and this will be longer than the time 
required to resolve frequencies to 1 Hz. 
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Another consideration in the design of active ranging 
detection systems is the requirement to accomodate doppler 
in the returning echoes. The input filter must have a 
bandwidth wide enough to pass echoes that are doppler 
shifted by the maximum amount anticipated during operation 
of the system, or the echoes will not be detected. If 
a filter of sufficient bandwidth to ensure this does not 
give a sufficient S/N ratio improvement, one may alter- 
natively use a bank of narrow band filters centered at 
various frequencies corresponding to different doppler 
shifts. This is illustrated in the diagram on the next 
page. 
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V. RANDOM SIGNALS , POWER SPECTRAL DENSITY, AND NOISE 



A. RANDOM SIGNALS 

Heretofore the signals discussed have been explicitly 
described "deterministic signals", such as a square wave, 
sinusoidal wave, sawtooth, etc. In describing a signal ex- 
plicitly as sume function of time, it is assumed that its 
amplitude and phase are known exactly for all time; i.e., 
past, present, and future. No real world signal meets these 
criteria, for then the signal would convey no new informa- 
tion, as it could be predicted exactly for any future time. 
Therefore, consideration must be given to random signals. 
Random signals of interest to the Antisubmarine VJarfare spe- 
cialist emanate, for example, from various machinery and flow 
noises. Noise in the ocean due to waves, wind, biological 
sources, seismic disturbances, etc., is a special case of 

random signal, and it is always present to some degree wheth- 

/ 

er or not a "target" is present. Processing equipment must 
be able to distinguish between the random signal of interest, 
unwanted signals, and noise. 

A random signal, unlike a deterministic signal, may only 
be described by its statistical characteristics. Its past 
history and its predicted future will be stated in terms of 
some average value and a given probabi lity that the signal 
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will be within certain limits at a specified time. This 
"given probability" is represented by a probability density 
function (pdf). 

What is a probability density function? A function is a 
descriptor, i.e. , it usually describes something. Density is 
a measure of quantity contained within a specified space. 
Probability refers to the chance, or likelihood , of the occur- 
rence of an event. Therefore, a probability density function 
describes the likelihood of some quantity being contained 
within a specified space. 

How does the probability density function relate to a 
random signal? Assume a random signal, x(t), which has some 
randomly varying voltage amplitude over time. 
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Plot the signal with a vertical scale representing voltage 
amplitude and the horizontal scale representing time. 




Suppose now that 1-volt represents any value between 0.5 and 
1.49; 2-volt represents 1.5 to 2.49; 3-volt represents 2.5 to 
3.49, etc. From the random, signal, a plot is now made with 
the vertical scale still representing voltage amplitude but 
v/ith the horizontal scale representing the number of times 
the signal is at that amplitude, or within that amplitude 
range. 

During time interval 1, the signal was present only in 
the 3- volt range and that one block is shaded in on the plot. 




c- rn f U i rc>- C 

During time interval 2, the signal is present in the 1, 



2, and 3-volt ranges. Those additional blocks are shovm in 



the next diagram. 
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During time interval 3, the signal is present in tne I 
and 2- volt ranges. Those additional blocks are shaded in. 
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The development should start to become apparent. Look 
now at time intervals 4, 5, and 6. The signal is present in 
the 3-volt range once, the 2-volt range twice, and in the 1- 
volt range three times. Shade those in. 
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To finish the plot, look at time intervals 7,8,9 and 10. 
The signal appears in the 4-volt range four times; in the 3- 
volt range three times; in the 2-volt range four times; in 
the I-volt range four times; and the 0-volt range once. 

Shade them in . 



The plot above represents the amplitude ranges that the 
signal v;as in for the entire time it was observed. The total 
number of shaded blocks represents the number of times the 
signal is within the various amplitude zones during the var- 
ious time intervals for the entire duration of the signal: 
the 29 blocks account for 100% of signal time. (Count them.') 
Each shaded block is regarded as a sample value. If the num- 
ber of sample values is denoted by n, then n = 29 in this 
example. The percentage of time the signal is in a 
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particular voltage range can be approximated by dividing the 



number of signal sample values 
total number of sample values, 

0-volt range: 
L-volt range: 

2- volt range: 

3- volt range: 

4- volt range: 



in that voltage range by the 

n. 

1/29 = 3% 

9/29 = 31% 

8/29 = 28% 

7/29 = 24% 

4/29 = 14% 



Total 100% 

The percentage of time the signal spends at other ranges 
can also be calculated. Suppose, for example, that in an ASW 
system there is a detector which detects only signals with 
amplitudes of 3 volts or more (the detector has a "detection 
threshold" of 3 volts). To determine the percentage of time 
the above random signal is detectable, it is necessary to de- 
termine the percentage of time it remains in the range of 3 
or more volts. The 3-volt range covers from 2.5 to 3.49 
volts. Take one-half the time the signal is in the 3-volt 
range and add it to the entire time the signal is in the 4- 
volt range to approximate the total time the signal is at or 
above 3 volts. 



(1/2) (7) + 4 = 7.5 Divide by n to get percentage. 
7.5/29 = 26% Therefore, the signal is detectable 

26% of the time. 

The accuracy in the determination of the percentage of 
time the signal is within any specified frequency range can 
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be increased by making the voltage amplitude ranges smaller 
and shortening the duration of the time intervals. The 
next diagram shows sample values for the same signal, but 
with amplitude ranges of 1/2-volt and the time interval dura- 
tion halved. It is constructed in the same amnner as the 
previous plot. A check to verify its accuracy will help you 
evaluate your understanding of the procedure. 




As amplitude range and the time interval length are made 
smaller and smaller, the plot approaches a continuous curve. 
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In most practical applications a formula for the continuous 
curve can be determined, and the corresponding function inte- 
grated over the amplitude range of interest to find the per- 
centage of time the signal was within that range. 

The probability density function of a signal represents 
the "probability" that a signal will be within an amplitude 
range. As the plot constructed above represents 100% of the 
time of the signal, so too does the probability density func- 
tion represent the total probability of the signal. It is 
scaled from 0 to 1, i.e., the area under the probability den- 
sity function curve is equal to 1. If the detector in the 
example just mentioned detects, with certainty, any voltage 
amplitude of 3 volts or more, and the signal spends 26% of 
the time at 3 volts or more, then the probability that the 
signal will be detected in the 3-volt or more range is 0.26. 
The probability density is expressed by a function which may 
be integrated to find the probability that the signal is in 
any prescribed amplitude range. For probability density func- 
tions used in many applied problems, the integration process 
can be avoided, sincd the results are readily available in 
standard tables. 

The probability density function (pdf) most commonly used 
for electrical signals is the "gaussian" or "normal" pdf. It 
is the familiar bell-shaped curve described by the equation 



- (x-m) 
2o^- 



p (x) = 



a/ 2tt 
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In the equation for p (x) , 



m = mean, or average, value, 

X = value of the signal, x(t), at time t, ^ 

a = standard deviation (a is the Greek lower case sigma) . 
The standard deviation, a, is a measure of the spread of the 
signal abour its average value, m. In the next figure we 
we show two gaussian curves with the same mean, but different 
values of o. The area under each curve is equal to 1.0. 





m 



For the gaussian pdf, about 
lies between m - a and m + a 
that the random signal, x{t) 
and m + lo is 0.68. 

m - la, P (x) = 0.68 

(Shaded area =68% 
of total area) 

mi 2a, P(x) = 0.95 
(Shaded area = 9 5% 
of total area) 



8% of the area under the curve 
which means that the probability 
assumes a value between m -la 
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m - 3a, P(x) = 0.997 
(Shaded area = 99.7% 
of total area) 




In the illustrations on the previous page, only 0.997 of 
the total probability has been accounted for. This is because 
there are no limits on a gaussian pdf; the range is unbounded. 
A signal may conceivably exist at any finite amplitude above 
or below the average value. The probability that the signal 
will occur outside the m - 3a range is only 0.00 3 (]_.0 - 0.99 7 
= 0.003). Limiting circuits are normally used early in the 
processing stage to prevent the rare occurrence of a large 
amplitude signal from damaging equipment. The important part 
of the signal, 0.997 of it, is still retained. 

The gaussian pdf is centered around an average value. 

For an electrical signal this is the average value of the 
amplitude of the signal. Look at the signal and its pdf 
side by side with the pdf oriented to the signal amplitude. 
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The value chosen for m is the statistical average of the 
29 sample values above. Each of the sample values was 0, 1, 
2, 3, or 4. The average, m, is equal to the sum of the indi- 
vidual sample values divided by the sample size, n. 

0+l+l+]+l+l+]+]+I+l+2+2+2+2+2+2+2+2+3+. . .+4. . . 

29 



or, since many values are the same, 
m 



_ 1(0) +9(1) + 8(2) + 7(3) + 4(4) _ . 

29 



_ 1 



n 



n 



Algebraically expressed, this is m = ^ where the 

means the sum of all the samples, x^ , from j = 1 to j = n, or 

x, + x„ + X., +. . . + X . If the number of sample values becomes 

very large, the sum of the products of the individual sample 

values and their probability of occurrence becomes a better 

choice for m. Increasing the number of sample values gives a 

more accurate description of the signal. If ra is calculated 

for the plot with the amplitude range and time duration halved, 

n = 74 and m = 2.095. The algebraic expression, if n is very 

large, is approximately “ 

m = (Xj )• P (Xj ) 



It has been shown how statistical characteristics may de- 
scribe a signal. The next problem is how to obtain signal 
characteristics that are statistically valid to identify a 
signal. 

The validity of statistical methods in determining the 
outcome of a process depends on repition of the process over 
and over again. With a random signal, x(t) , this repetition 
may be accomplished by using an "ensemble" of random signal 
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generators with each signal generated having the same broad 
statistical characteristics. If the individual outputs from 
this "ensemble" are measured simultaneous Iv at some time, t, 
and averaged, the result would be an "ensemble average", in- 
dicated by x(t). This "ensemble average" is identical to, 
and has all the properties of, a statistical average. How- 
ever, most signals of interest in an ASW environment originate 
from a single source. Therefore, only the time average of 
the signal, ^(t)^, may be obtained. 

' The time average, ^ (t^/ and the ensemble average, x(t), 
are not always the same. For example, suppose the statisti- 
cal characteristics of the signal from the random signal gen- 
erators in the ensemble vary with time. Such a variation 
would not be reflected in measurements made at a fixed time, 
tj^. The ensemble average at time tj^ would be different from 
the ensemble average at time t 2 * A random signal with these 
characteristics is called "non-stationary " . If a detector is 
used to detect the random signal in our example above, what 
is the effect on the detector if the random signal is "non- 
stationary"? The average, m, was calculated as 2.14. Suppose 
that m varies with time, and that over another time interval 
m is calculated to be 1.14. If the signal over the second 
time has the same characteristics as over the first, the 
effect on detection would be the same as raising the thresh- 
old to 4 volts over the first interval, as indicated in the 
next diagram. 
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The probability that x = 3 volts or more would become 2/29 = 
0.069 vice 0.26 for m = 2.14! If the average is unchanged 
and the standard deviation increased by a factor of three, 
the detector would be observing only 0.68 of the signal as 
opposed to 0.997. Therefore, if the random signal is not 
"stationary" (i.e., its statistical characteristics do not 
vary with time) or at least "stationary" for the processing 
time of the detector, the detector will be very inefficient. 
Even though the signal may be stationary, the time and ensem- 
ble averages may still be unequal. 

Fortunately, many random signals that mjst be processed 
in ASW have time and ensemble averages which are identical, 
perhaps not for all time, but at least for the processing 
time required. Signals with their time and ensemble averages 
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equal , <^x ( t = x(t), are said to come from an "ergodic pro- 



cess", or to have the property of "ergodicity " . An "ergodic 
process" is also stationary. Under the assumption that a 
signal from an ergodic process means the time and ensemble 
averages are equal, the signal is stationary and the time 
average will have all the properties of a statistical average. 

Having extablished that a time average will be statisti- 
cally valid if it is from an ergodic process, how do we ob- 
tain the other statistical characteristics which describe a 
signal? All acoustic detectors convert the signal (and, un- 
fortunately, the noise, which v;ill be addressed later) into 
some time- varying voltage or current, x(t). 



city, is equal to the average value, m. This average value 
is both the value about which the signal's gaussian pdf is 
centered and, as the time average, represents the "dc" compo- 
nent of the signal (direct current, also called the steady- 
state or the non- time-varying component). If m = 0, there 
is no "dc" component. If the signal, x(t), is impressed a- 
cross a one-ohm resistor, the instantaneous power dissipated 
is x^(t). The time average, ^^(t)^, becomes the "second 
moment", which represents the total average power, x^ , of the 




X (t) 



The time average of the signal. 




assuming ergodi- 
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electrical signal. The "second moment" differs from the 
"average value squared", m^ , which is the dc power, or the 
power in the non-time-varying component: m^ is obtained by 
squaring after averaging, whereas is obtained by squaring 
before averaging. 

As an illustraion of these values, consider a random 
signal , x (t) . 
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The following drawings illustrate computation of and 
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There sre two other measures which identify a specific 
gaussian pdf besides the average value, m, about which it is 
centered. A random signal is composed of a dc component and 
an "ac" component (alternating current, or the time varying 
component of the signal). The total power, , must be the 
sum of the dc power, m^ , and the ac power, a^. Though repre- 
senting the ac power of the signal, is the "variance" of 
the gaussian pdf. The variance is equal to the second moment 
minus the average squared, = x^ - rn^ . For an electrical 
signal the variance is the total average power minus the dc 
power. As the standard deviation, a, represents the square 
root of the ac power, it is more commonly known as the root- 
mean square (rms) value. We summarize these definitions 
briefly. 

Given an ergodic random signal, x(t); 

1. the average (mean) value, m, is its dc component; 
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the average-squared, , is its dc power; 

3. the second moment, , is its total average power; 

4. the variance, - m^, is its ac power; 



5. the standard deviation, a, is its rms value. 

In "Fourier Transform Properties", the concept of correla- 
tion of a signal was discussed. It is considered again here 
because of its unique properties, and because of its relation- 
ship to the "power spectral density" of a random signal. 

Those who feel they have not grasped the concept of correla- 
tion should review that section before proceeding. 

The autocorrelation function, R(t), was defined previously 
as the time average, R(t) = ^x(t)»‘x(t + For an ergodic 

random signal the autocorrelation function becomes: 



The value of R(t) is a function of the time delay, t. At x=0 , 
R(t) is a manimum: 



As increases, R(t) decreases, until at t= infinity, R(t) 
is a minimum: 

R(<») = m^ = dc power; 

R(0) - R('») = = x^ - m^ = ac power. 

The autocorrelation function is also a measure of the 
"time coherence" of the random signal. If t is small, x(t) 
and x(t+x) will be close together in time. Therefore, the 
presence of one will have some effect on the presence of tne 
other, a condition known as "statistical dependence". If the 




R(0) = x(t) • x(t+0) = x^ = total average power. 
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presence of x(t) has no effect on x(t+T), such as happens if 
T is large and x(t) is random and non-periodic, then x(t) and 
x(t+T) are said to be "statistically independent". If for a 
signal with a gaussian pdf the absolute value of the auto- 
correlation function minus the dc power equals zero, |R(T)j - 
= 0 , then x(t) and x(t+T) are statistically independent 
and "uncorrelated". Statistical independence simplifies the 
calculation of the autocorrelation function. Whether or not 
two signals are "correlated" becomes important in the calcula- 
tion of their average powers. 

Assume that a signal has been formed by the addition of 
two separate signals: z(t) = x(t) + y(t). The correlation 

function of z (t) will be of the form 

(x+y) (x+y) = x^ + xy + yx + y^ , and is equal to: 

Rz (t ) = ^x (t) ‘ X ( t+T ^ ^ (t) • y (t+T )^ + ^ ( t) • y ( t+T )^ + ^ ( t) • x (t+T 

^ + \x'"> 

where the last two terms are the "cross correlation temis" of 

x(t) with y (t) and y(t) with x(t). As noted previously, the 

above correlation functions represent the average power in 

the signals, and can be written as: 

p=p + P + P + P. 
z XX yy xy yx 

If the two signals, x(t) and y(t), are correlated over the 
time delay, T, then the total signal power, S, is the sum of 
the signal power in signal x(t), the signal power in sig- 

nal v(t), S , and the cross correlation of their combined 

yy 

signal powers, and S : 

'^y y^ 
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s.-=s +S +s +s 

z XX yy xy yx 

If the signal is the sum of mutually un cor re la ted signals, 
as in the case of noise un correlated over the time delay, t, 
then the total noise power, N, is the sum of only the noise 
power and The cross correlations are equal to zero: 

vCX vy “ — — — — . — — — 

N = N + N ; (N = 0 ; N = 0 ) . 

2 XX yy' ' xy ' yx 

For a signal which is the sum of mutually uncorrelated 
signals, the correlation function is the sum of only the auto- 
correlations (the cross correlations equal zero) , and the av- 
erage power is the sum of only the average powers of the 
individual signals (the cross correlated powers equal zero) . 



B. POWER SPECTRAL DENSITY 

Til us fai. Oil ly the time domain aspect of random Signals has 
been addressed. The frequency domain is represented by the 
power spectral density, which was introduced in "Filters and 
Bandwidth Requirements". The device used to transfer between 
the time and frequency . domains is the Fourier Transform. The 
Weiner-Kin chine Theorem states that the power spectral density, 
G(f), and the autocorrelation function, R(t), are Fourier 
Transforms of each other: 



'^oo 

G(f) = ^ [R(t)] = r x(t) x(t+x)e ^^^dT . 

And, from the Duality Theorem, 

if R(x) -^-^G(f), then G(x)-^-^ R(f). 

As can be seen above, the random signal may be described 
in either the time domain," R(x) , or the frequency domain, G(f) , 
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with free interchange between the two. 

Since the power spectral density is the Fourier Transform 
of the autocorrelation function, the Fourier Transform Theo- 
rems may be applied to the power spectral density. The sys- 
tems designer uses the properties expressed by these theorems 
to design an efficient detector/signal processing system. The 
ASW specialist who has knowledge of these properties can bet- 
ter understand how his system is designed to operate and 
therefore, utilize it more effectively. The duality theorem 
was already utilized. The other major theorems applicable to 
power spectral density will be explained along with examples 
of their usage. 

1. Frequency Translation Theorem 

If a signal, x(t), is bandlimited in W < f^, and a 

new signal, y(t), is formed by multiplying x(t) by 

cos (to t+(|)) , y(t) = x(t) cos (to t+(t>) , then the correlation 
c c 

function of y(t) is equal to one-half the correlation 
function of x(t) costo t, that is R (t) = (]./2)R (T)costo t. 

C y K 

In addition, the power spectral density, G(f), is 
Gy(f) = (l/4)G^(f-f^) + (l/4)G^(f+f^) . 

Multiplication in the time domain by costo^t becomes transla- 
tion in the frequency domain, as the power spectrum is shifted 
up and down in frequency by The condition that the band- 

width, W, be less than f^ is necessary to ensure that G^(f-f^) 

and G (f+f ) do not overlap, otherwise aliasing, as described 
X c 

in "Filters and Bandwidth Requirements'; would result. Notice 
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also that in multiplication by cos (o)^t+(|)) , the phase factor, 
<J), is lost in translation to the power spectral density. 
Therefore, the original random waveform may not be recon- 
structed from knowledge of the power spectral density. How- 
ever, by returning to the time domain via the inverse Fourier 
Transform, the autocorrelation function, R(x), with its sig- 
nificant statistical properties, may be found. 

Frequency translation is very useful in signal processing 
Suppose that due to space, weight, or power limitations it is 
impractical to locate the detector and the processor together 
The low frequency (0 to ]000 Hz) sound may be detected, trans 
lated to radio frequency (RF) (HF =2-30 MHz, VHF = 100 - 
200 MHz, UHF = 200 - 400 MHz) and broadcast to the processor. 
In the same manner, one processor can handle many detectors 
if each detector's signal is translated to a separate fre- 
quency band within the RF range (Detector 1, frequency f to 
f 2 ; Detector 2, frequency f^ to f^; etc., all within the UHF 
band) . The problem could be one of processing signals of 
many different frequencies. Rather than having many differ- 
ent processors, each processing only one signal, it may be 
better to have one processor operating at a fixed frequency. 
All of the signals could be translated to the processor fre- 
quency for processing. If the processor is very elaborate 
and/or very expensive, this method becomes desirable. 

2 . Integration Theorem 

1 

If y(t) = J x(t')dt', then G , (f ) = — G (f ) . 

-too y (2Trf)^ ^ 
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The multiplication factor, that the low 

frequency components of the power spectral density will be 
enhanced by integration. The ocean acts as an "integrating 
filter", in that it allows the low frequencies to propagate, 
while rapidly attenuating the higher frequencies. This is one 
reason why, in an attempt to get long detection ranges, active 
sonars and passive processing systems have increasingly ex- 
ploited the lower frequencies. 

3. Differentiation Theorem 

This theorem is offered without an example since it 
follows so naturally the integration theorem. 

If Y(t) = d[x(t)]/dt, then G (f) = ( 27 Tf ) (f) . 

Y ^ 

The multiplication factor, (2Trf)^, indicates that the high 
frequency components of the power spectral density will be 
enhanced by differentiation. 

C. NOISE 

At the beginning of this section, noise was mentioned as 
a "special case of random signal which is always present to 
some degree whether or not a 'target' is present". In com- 
munications, noise is often defined as any electrical inter- 
ference, or unwan ted signal . Signal then, in addition to the 
characteristics previously used to describe it, has the qual- 
ity of being wanted, sought after, emanating from a source 
of interest, and containing intelligence. In ASW, signal is 
always associated with "target". Noise, on the other hand. 
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is anything which interferes with, or tends to. mask, the 
signal. Often "one man's signal is another man's noise" as 
any patrol aviation Tactical Coordinator whose buoy field has 
just been penetrated by a "pinging" Destroyer will verify. 

An acoustic detector converts both signal and noise into a 
time-varying electrical voltage or current. An examination 
of some of the characteristice of noise will help in under- 
standing how systems may be designed to reduce it. 

Another way to define noise is simply as that quantity 
observed in the absence of signal. Noise may come from a var- 
iety of manmade or naturally occurring sources. Systems are 
designed to eliminate as much noise as possible, but some 
noise will inevitably remain. The thermal motion of electrons 
in the conducting media of the com.ponents of a sysrem, for 
instance, is one unavoidable cause of electron noise. Ther- 
mal noise is interesting in that its power spectrum is con- 
stant over a wide range of frequencies. It is designated 
'‘^white noise" by analogy to white light, as all frequency com- 
ponents are present in equal amplitudes. The amplitude of 
thermal noise has been proven to have a "gaussian" pdf. It is 
therefore referred to as "gaussian white noise" ("gaussian" 
amplitude distribution; "white" frequency distribution). If 
"gaussian white noise" is played over a loudspeaker, it sounds 
dull and monotonous, somewhat like a waterfall. The subtle- 
ties of its random variations are hidden from the human ob- 
server. 

The power spectrum of gaussian white noise, with zero 
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mean, is G(f.) ^ » where n is pov;er density in watts per 

Hertz. By Fourier transformation, the autocorrelation func- 



tion is 



e^“^ df - (^)6 (t) , where 6 represents 
the "delta function" introduced previously. 
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It is apparent from the R(t) plot that for any time delay, t, 
other than zero, the autocorrelation function is zero; there- 
fore, any two different samples of a gaussian white noise sig- 
nal are uncorrelated and statistically independent. 

The signals of interest in ASW are frequency limited, 
i.e., there is some finite frequency band within which the 
signal exists. Filtering in order to improve detection of 
that signal was discussed earlier. If gaussian white noise 
is filtered, its frequency components are naturally those of 
the filter bandwidth, but the amplitude distribution remasins 
gaussian. The output power spectrum of white gaussian noise, 
filtered by an ideal filter, is a rectangular function. 



^y ~ ^ ’ where II indicates a rec- 
tangular function and B represents the filter bandpass. The 
autocorrelation function is the inverse Fourier Transform of 
the power spectral density. As indicated in "Fourier Trans- 
form Properties", the Fourier Transform of a rectangular 
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The above figure shows that filtering has produced the follow- 
ing results. 

1. The power spectrum, though no longer white, is con- 
stant over the finite frequency range of the filter. 

2. The output power is finite, ~ ^ (2B) = riB. It 
varies linearly with bandwidth, B. 

3. The output signal is correlated over time intervals 
of about 2 ^* 

How is this knowledge used in signal processing? 

Item (1). The noise in which a signal of interest must 
be detected is often not gaussian white noise. However, 
through judicious choice of a bandwidth, its power spectral 
density may be constant over that range. Otherwise, a pre- 
whitening filter may be used to make the noise power spectral 
density constant over the bandwidth. Therefore, after filter- 
ing it may be treated as gaussian white noise. 
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Item (2) . The noise power after filtering is = r|B. 

It is a direct, linear function of bandv;idth, B. By narrow- 
ing the bandwidth, the amount of noise power present is re- 
duced. Therefore, the signal power is enhanced in relation to 
the noise power. This is crucial in a power, or energy 
detector. 

Item (3). The bandwidth of the filter determines the min- 
imum time delay, t, in a correlation detector. It must be 



will be correlated. If the noise is uncorrelated, which it 



lation detector, thus decreasing in relation to the signal 
power. 

The bandwidth of the filter also determines the finite 
integration time to be used in determining the time average. 



A review of the development of bandwidth versus integration 
time in the "Filters and Bandwidth Requirements" section may 
be helpful in the understanding of these relationships. 



of random signals, power spectral density, and noise. The 
intent has been to provide a basis for understanding their 




v;ill be for t> only its average power will add in a corre- 



The time average is indicated by (brackets) and is defined 



as : 




,T 



Replacing T by ^ and by we obtain 




0 



This section has been an attempt to describe the concepts 
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individual characteristics. A grasp' of these characteris- 
tics is necessary to understand the methods of signal process- 
ing. The most common methods in current use, DELTIC, energy 
detectors, correlation detectors and beam formers, will be 
presented in following sections. 
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VI. THE DELTIC 



In the section on filters and bandwidth, the relation- 
ship between integration time and frequency resolution 
was discussed. It was shown that one must integrate for 
longer periods of time as the frequency resolution 
requirements become stricter. In many cases, however, 
the integration time required to give the necessary 
resolution is much too long to be of tactical use. In 
order to circumvent this, several processing schemes 
have been developed which in effect give more resolution 
for a given integration time. One of the more common 
in use with military processors is known as the DELTIC 
( DEL ay Time Compressor) . Let us investigate the manner 
in which this processor accomplishes this feat. 

It has been shown that processing the signal as it 
is received by the detector will not give any greater 
resolution. The DELTIC takes the signal as received, 
and records it on a magnetic tape, disc, or drum. It 
does this in a special sequence, however. To see this, 
let us take a simple signal and observe how it is processed. 




For simplicity in showing the operation, an 8-sainple 
point "time window" will be used. In actual practice, 
1000 or more sample points are commonly used. I£ the 
time between sample points is T, the time window observed 
will be 8T = T. Therefore if the DELTIC is not used, 
the best resolution we can get is 1/8T, and it would 
require Sx seconds to get this resolution. The physical 
layout of a DELTIC is shown in the next diagram. 




The drum is rotated at a rate such that it makes one 
revolution plus one sample space during the time period 't. 
In this case, with 8 sample points it makes 1.125 rev- 
olutions per 't' seconds. The read-out head continuously 
reads out the sample values as the drum rotates under it. 
The input gate injects a new sample point once each 
1.125 revolutions in such a way that the sample point 
which has been on the drum the longest time is replaced. 
Starting with a blank drum, let us examine the process 
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When the drum has made 1,125 revolutions [in sec.) 
the input gate inserts sample point #2. 




After 1.125 more revolutions, the 3rd sample point 
is inserted. 
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The remaining spots are filled similarly until all 
8 are fill'ed. The output during one revolution of the 
drum then is : 
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When the next sample point is inserted after 1.125 revo- 
lutions, it displaces the first. 
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Once the drum is filled (which takes 8T sec.), a new 
value replaces the oldest one every T seconds. If we 
look at the output from the read-out head, we see that 
there are 9 values read out each T seconds (8 values/ 
revolution, 1.125 revolutions/r seconds). 
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The output of the DELTIC each T seconds is thus the 
same as the inputs over the last 8^ seconds. The DELTIC 
operation described above can be illustrated by the 
following analogy. If we record 80 seconds of conversation 
on tape, and then speed up the playback so that the 
conversation takes only 10 seconds, we have compressed 
the conversation just as a DELTIC would. The same 
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information is in the 10 second playback as was in the 
80 second conversation, but it only takes l/8th the time 
to hear it. We also notice that this compressed conver- 
sation has been shifted in frequency by the compression, 
and normal voices sound like Donald Duck. This illustrates 
the practical effects of the Fourier scale change theorem. 

By using the DELTIC ^ve then are able to observe the 
preceding 8T seconds* worth of signal during each interval 
of length T seconds. The frequency domain plots of the 
original signal and the output of the DELTIC show that 
the original spectrum has been shifted in frequency and 
stretched out. 
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The Fourier scale change theorem shows that we have 
multiplied the time domain characteristics o£ the signal 
by the factor 1/9, and so the frequency domain character- 
istics of the signal have been divided by this factor, 
thus effectively multiplying the frequency plot by 9. 

The "compression factor" is the ratio of the original 
sample time, T, to the compressed sample time, T'. 

Thus in this case it is the ratio T to T/9, or 9. Note 
that not only is the center frequency and bandwidth 
multiplied by this factor, but also the minimum resolvable 
frequency difference. Without compression, we would 
integrate for T seconds to obtain a resolvable difference 
of 1/T. Upon compression, this becomes 9/T. At first 
glance, this does not seem to i)elp, since the resolution 
has not been improved. Why use the DELTIC if it doesn't 
help? 

The benefit in using the DELTIC is that we can obtain 
this resolution not by integrating over T seconds, but 
over T/9 seconds once the DELTIC has been filled. Many 
spectrum analyzers "scan" the frequency band of interest 
by using a single band-pass filter of bandwidth equal 
to the minimum resolvable frequency difference of the 
processor, and stepping its center frequency across the 
band one step each integration period. Thus, if the 
filter bandwidth is 1/9 the bandwidth of interest, without 
the DELTIC it would require 9 integration periods (T) for 
one scan of the band. With the DELTIC, each integration 
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period is T/9, and thus we can scan the entire band with 
the same resolution in what was previously one integration 
time, T. This is quite an advantage for a band of 
several hundred hertz and a required resolution of 0.5 Hz. 

The DELTIC need not be designed for insertion of a 
new sample value every revolution, but only a new value 
approximately every 10 revolutions. In this way one can 
compress the signal even further by having the drum 
revolve more than once every q* seconds. 

Instead of a magnetic drum, disc, or tape memory, 
a digital shift register may be used. The shift register 
is stepped past the read-out gate in such a way as to 
read the entire register each x seconds. The values 
are then recirculated. Each time a complete circuit is 
made, a new value is input to replace the oldest, and 
the process continues. 
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SoiiiC' systems which use tlie DOLTfC in processing ;ire: 

a. . AN7AQA-7 

b. A\7BQR-13 

c. A.\/SQS-26 

d. AN7BQR-20 

The illustrations on the following pages show a side by 
side comparison of an input signal and the resultant DF.LTIC 
output. Notice tlie gradual build up of tlie CELTIC output as 
time intervals of the signal are introduced. .Note also that 
as time intervals continue to be introduced, the CELTIC drops 
out the oldest information and represents only the most re- 
cent time intervals. 



110 



A 




As the signal is "stepped” into the DELTIC, the output 
of the DELTIC represents the compressed signal. Plots A to 
F indicate the buildup of the DELTIC. Plots G to J show the 
oldest information dropping out and only the most recent in- 
formation being represented. Plot K shows the complete time 
history of the signal and its DELTIC output. 
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VII. ENERGY DETECTION 



In the discussion of random signals it was noted 
that in order to apply statistical analysis to signal 
processing, the signals must be "ergodic”. One property 
of ergodic signals is statistical stationarity , that is, 
the statistical properties of the s ignal are constant 
in time. In reality, the signals that must be processed 
are not strictly stationary. The sea noise varies with 
the time of day, changes in weather, etc., but the var- 
iation is over a long time period. This fact enables 
us to apply the statistical analysis approach to signal 



that the variation in signal properties is negligible. 

Let us consider how to apply statistical methods to 

detection problems. For simplicity, take the case of 

Gaussian noise, where the probability distributions of 

amplitude for noise alone and noise plus signal are 

2 

Gaussian Avith equal variance, (S , but different means 




(average values) 
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If we now have a processor which compares the amplitude 



of the input to a "threshold" level set into it, 
indicates "signal present" if the input is above 
"signal absent" if the input is below .threshold, 
we can apply statistical theory to determine the 
ating characteristics of the processor. 







and 

threshold , 

then 

oper- 



The area under the S + N curve above the threshold 
is proportional to the probability that the amplitude of 
the input will be above the threshold when the signal 
is present. Thus, this area is proportional to the 
signal detection probability, P(d). 







The area under the N curve and above the threshold 
is proportional to the probability that the amplitude o£ 
the input will be above the threshold, even when the 
signal is absent, because o£ the background noise. Thus, 
this area is proportional to the probability o£ £alse 
alarm, P(£a) . . . , 



a 



I£ the threshold is varied, the probabilities change. 
I£ the threshold is raised, the probability o£ detection 
decreases, but so does the probability o£ £alse alarm. 
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If the threshold is lowered, P(d) increases, but 



so does P ( fa) . 




The ratio of P(d) to P(fa) does not remain constant, 
however. Plotting P(d) versus P(fa) as the threshold 
varies, we obtain the following resulL. 




This is known as a "Receiver Operating Characteristics" 
curve, or ROC curve. Using this curve we can determine 
where to set the receiver threshold to optimize P(d) 
and P(fa) 
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As the threshold is reised, P(d) .does not change 
much in comparison to the rate at which P(fa) is 
decreasing, until the "knee" of the curve (between point 2 
and point 3) is reached. 





As we continue to raise the threshold above this 
value, P(d) starts to fall off rapidly, while P(fa) 
changes very little. By deciding what maximum P(fa) 
can be tolerated, we can set the threshold to yield the 
maximum P(d) possible for that particular P(fa) by 
referring to the ROC curve. Conversely, if a minimum 
P(d) is required, the ROC curve indicates the minimum 
•P(fa) that must be tolerated. 

Until now the effect of changing the mean of the 
S + N distribution with respect to the N distribution 
would have on P(d) and P(fa) has not been discussed. 

If the mean of the S + N distribution is raised, then 

the ratio P(d)/P(fa) is larger for all threshold settings. 
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The result is shown on the ROC curve in the following 
diagram . 




We define the "detection index" as the quantity 

/C , which is a measure of the difference 
between the distributions. A large d value indicates a 
large separation between the distributions, and thus a 
curve with more of a "knee". 
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P(d) 




P(£a) 



This diagram shows that if the presence of a signal 
drastically changes the input distribution mean, then 
we can reduce the P(fa) significantly while still 



does not drastically change the mean, large values for 
P(d) are possible only at the expense of a large P(fa) '• 

We hav'^e discussed thresholds, and the way the selection 
of threshold values affects the probability of detection 
and the probability of false alarm, but have not discussed 
how these thresholds are set. One point that may not 
be readily apparent is that in detection systems in which 
the human operator is involved in the decision as to 
presence or absence of a signal, his perception of the 
display must be included in computation of the threshold 
value. Most systems in use operationally involve an 
operator in the decision process. The SQS-23 sonar is 
an example. The operator can vary the threshold to an 
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extent by adjusting the brightness and intensity of the 
display and the gain control, but a large consideration 
in determining whether or not a target is present is his 
ability to distinguish betw^een target and noise as pre- 
sented on the PPI scope. 

The operator reading grams on an AQA-7 is another 
example. He can control the gain, or the grey value of 
the background on the paper, but his perception of the 
distinction between target and noise determines the 
threshold. For this reason it is difficult to determine 
the ROC curves for systems with human operators, because 
it is difficult to determine the actual threshold setting. 

The threshold will vary from day to day, or even hour 
to hour, vn th the same operator and the same hardwaie 
settings because of the operator's state of mind, physical 
condition, fatigue, etc. This is what makes quantification 
of the detection characteristics for a system difficult. 
Completely automatic systems are much easier to evaluate, 
as one can determine exactly what the settings are. 

In the preceding discussion of energy detection and 
thresholding we have seen how a system's detection capability 
can be determined from knowledge of the characteristic 
distributions of noise and signal plus noise. Implicit 
in this treatment is the assumption that the distributions 
of noise and signal plus noise are exactly known. If these 
distributions are exactly known we can specify the per- 
formance of the system for any given threshold setting 



120 



in terms of i*(d) and P(fa) . If we do not know these 
distributions exactly, our analysis will not exactly 
describe the system performance. Thus, the accuracy of 
our prediction of system performance depends on the 
accuracy with which we can describe these distributions. 

In designing real world systems this is one of the problems 
facing the designer , since the distributions of noise and 
signal plus noise are not constant in time (over any 
lengthy period) or with geographical position. For this 
reason many systems do not always perform as well as 
the specifications indicate. Only if the ambient conditions 
are the same as those for which the system was designed 
will it perform as predicted. 
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VIII. CORRELATION DETECTION 



In the section on random signals, power spectral 
density, and noise, it was shown that correlation can 
effectively improve the signal-to-noise ratio for cor- 
related signals due to the cross product terms in the 
calculation of the correlation function. How can this 
information be used to detect signals in noise? Let 
us study an active echo ranging system as an example. 

In this type of system a signal is transmitted, and 
then the returning signals are processed. The range 
of the reflecting object can be determined from the echo 
return time. In a noise- free environment, there are 
fev^ problems in this system, but in most real situations 
the ambient noise may "bury" some of the returns. By 
using correlation, we can improve the signal-to-noise 
ratio and recover these echoes. 

This system retains a replica of the transmitted 
signal and then inputs the received signal into the 
correlator. The output of the correlator is large when 
an echo is received, and small when only noise is present 
What does doppler do to the correlation of returning 
echoes? Let us look at a burst of CW signal transmitted, 
and the effect of doppler shifting on the correlation 
function. If the echo is not doppler shifted, the output 
of the correlator will be the autocorrelation of the 
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transmitted signal at a time , which is related to 
the range o£ the reflected object. 




transmitted 

signal 




returning 

echo 




output 



As the returning echo is shifted in frequency by 
doppler, the output of the correlator becomes the cross- 
correlation of two closely related signals, but its value 
is never as high as the autocorrelation. The cross- 
correlation amplitude and the amount of delay over ivhich 
it is correlated decrease. (In fact, the reduction results 
from tiiO decrease in the delay over w’nicii Liie signal is 
correlated with its doppler shifted echo.) 




transmitted doppler shifted correlator • 

signal returning echo output 



A plot of the correlator output for a number of different 
doppler shifts on the same diagram, that is, a *t, h£ , 
correlator amplitude plot, shows the effects of doppler. 
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For a given echo, this plot gives a solid surface 
relating Af, 't, and amplitude. 




If we look down at the T, Af-plane, and plot the 
contour line where the amplitude of the function is 
50 per cent of the maximum, we obtain the following 
resul t . 







This diagram is commonly known as an "ambiguity 
diagram", and the reason for this will become apparent 
from examination of FM pulses. First, however, let us 
consider a long CW pulse to see how doppler affects the 
correlator output. 

The section on correlation and convolution shows 
that the width of the correlation function is related 
to the duration of the time-domain signal. 



expected to respond over a long range of values of 
What about doppler and its effect on a long CW pulse 
being correlated? Observation of two signals at slightly 
different frequencies shows that the longer the pulses 
to be correlated, the less the correlation, since the 
signals will be almost correlated over only a few cycles, 
but get further and further out of phase as time progresses 



As the difference in frequency between the signals 
increases, the number of cycles over which the two are 
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almost correlated decreases rapidly. Thus for long 
CW pulses, the correlator response decreases rapidly 
with doppler shifts. The ambiguity diagram for a long 
CW pulse is shown in the following figure. 




The length L is proportional to the length of the 
CW pulse. is the delay corresponding to the actual 

lauge of the target. 

For a short CW pulse, the length of the correlation 
peak 50 per cent contour will be shorter, but the range 
of doppler shifts over which the signals remain essentially 
correlated is greater. 




Tlius the ambiguity diagram for the short CW pulse 
has the following appearance. 
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These two diagrams show why long pulse CIV is good 
for searching and doppler determination, and short pulse 
CW is good for range resolution. A sonar pulse of 120 
milliseconds will ensonify a region of water 180 meters 
long. In other words, the spatial length of the pulse 
as it travels through the water is ISO meters (1500 meters/ 
second velocity x 0.120 seconds = 180 meters, the distance 
the forward edge of the pulse has traveled during the time 
the rest of the pulse was being generated) . The cor- 
relation peak 50 per cent contour, as shown on the 
ambiguity diagram, will be proportional to the pulse 
length T, and for CW this turns out to be a one-to-one 
ratio. Thus, the correlation will be over a period of 
120 milliseconds, which means that the target range can 
be determined only to the nearest 180 meters. 

In the case of the short pulse of, say 5 milliseconds, 

/ 

the pulse spatial length is 7.5 meters, the correlation 
peaks over 5 milliseconds, and the range resolution is 
7 . 5 meters . 
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The effect of doppler on the amb.iguity diagram shows 
that the same type of development can be made to relate 
frequency resolution of the two different pulse lengths. 
Because the long pulse remains correlated over very 
limited doppler shifts, the frequency of the returning 
echo can be determined rather accurately by using a bank 
of correlators comparing the echo to replicas with varying 
amounts of doppler shifts. 




On the other hand, use of the same bank of correlators 
with a short pulse produces a response from each of these 
correlators, because the short pulse remains correlated 
over a large range of doppler shifts. 
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If we use a pulse in which the frequency is changed 
during transmission, which is commonly referred to as 
FM slide, several differences appear. 




Because of the frequency change during the pulse, 
the returning echo correlates with the transmitted replica 
over a short range of timOj^C, even for long pulses. 

Thus, the ambiguity diagram is very narrow in the X 
direction. On tlie otlier Imiid, if a lojig pulse i:? used, 
the Af range of correlation is narrow. One difference 
between FM and CW pulses becomes apparent here, however. 

In the case of FM pulses, a doppler shift will tend to 
shift the entire pulse, in frequency, but the slide will 
remain linear (provided it was originally linear), and 
thus it will still tend to correlate but at a different 
value of T. This behavior and its effect on the ambiguity 
diagram are shown in the next drawing. 



129 







“DoPPueR 

5.uiPTeo 

RSTURM 




Note that the FM slide pulse has some of the char- 
acteristics of both long and short pulse CW. At any 
given doppler shift, the correlation peak itself is 
narrow, as is the peak for short CW pulses. Thus, range 
resolution is good. 
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Note also that at a given T the range of doppler 
shift over which the signal is correlated is narrow, 
on the order of the range over which the long CW pulse 
is correlated. 







The reason this diagram is called an "ambiguity 
diagram" can now be seen. The correlator output, that is 
the T-amplitude plot, shows the following. 
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Was the peak at caused by a target exhibiting 
no doppler at range equivalent to or was it caused 

by a target with 5 knots o£ doppler at range - A? 
The output looks the same in either case. There is a 
"range ambiguity" present. 




One way of resolving this is to transmit two FM slide 
pulses, one with an "up-slide", (increasing frequency 
with time), and one with a "down-slide" (decreasing frequency 
with time). By then using two correlators, one with an 
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”up-slide" replica and one with a "down- si ide” replica, 
the ambiguity can be resolved. The next ambiguity diagram 
shows the output from both correlators for a given range 
target . 
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When compared on a T-amplitude plot, the correlator 
outputs for a target exhibiting the doppler shift shown 
on the ambiguity diagram by the dotted line appear as 
in the following diagram. 
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By noting which peak appears first, we can determine 
if the doppler is up or down. If the "up-slide" peak 
appears first, then the returning echo has been shifted 
down in frequency, and thus correlates at a delay, , 
earlier than it would if not doppler shifted. This 
indicates down doppler. 







ECHO 



If the "down-slide" peak appears first, the opposite 
is true, which indicates doppler. 
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Echo 
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The true range of the target is midway between the two 
correlation peaks, regardless of the amount of doppler 
shift present. If the target has no doppler, the two 
peaks coincide. 





The separation between the up-slide and down-slide 
correlation peaks is proportional to the amount of doppler 
shift in the returning echo. If a suitably calibrated 
display of the two outputs is available, the amount of 
doppler may be read directly from the presentation. 

Thus, by employing FM slide pulses and suitable 
processing, one can obtain the range resolution of short 
CW pulses and the doppler discrimination of long CW pulses 
from a pair of long FM pulses. 
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Another type of pulse which provides the "best of 
both worlds” with respect to doppler and range is known 
as the "pseudo- random noise pulse”, which consists of a 
"coded” pulse with amplitude or frequency varied in a 
pseudo-random fashion. Take, for example, the case where 
the amplitude (or frequency) of the signal is not a 
simple sinusoid, but a binary coded function. 




This pulse is correlated over a very narrow range 
of “t , giving good range resolution, and it is also cor- 
related over a very narrow range of doppler shift because 
of its pseudo- random character (Truly random signals 
correlate only at a single value of 'T.) The ambiguity 
diagram for a pseudo- random pulse is shown in the next 
figure . 




Note that because o£ the correlation over narrow 
ranges of A£ and T, no range ambiguity is present. The 
problem associated with this type of pulse, however, 
is tlia L tho preseiice of noise in the returning echo 
causes marked degradation in the correlation of the 
return, and thus degrades the performance of real 
systems. For this reason the pseudo- random noise pulse 
is useful in applications where the signal-to-noise 
ratio is high, such as close range tracking, bathymetry, 
etc. , but its usefulness is limited in long range search 
and tracking, where the signal-to-noise ratio is much 
lower . 

Up to this point we have covered the use of correlation 
in active systems, and how the use of FM slide pulses 
combined with the proper correlation processor enables 
us to determine the amount of doppler in an echo directly. 
Let us now examine some applications of correlation in the 
passive detection field. 
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One difference between active and passive correlation 
detection systems becomes apparent very quickly. That 
is, in the active system we compare the returning echo 
with a replica of the transmitted pulse. In other words, 
the signal we are seeking is a known quantity (a replica 
of the transmitted pulse). In a passive system, however, 
no pulse is transmitted, therefore we have no replica 
for comparison with the incoming signal. Therefore we 
cannot utilize correlation in the same manner; however 
it is still of use, particularly in passive direction- 
finding systems. The utility of passive direction- finding 
systems can be seen by examination of a typical system 
and its operation. Consider a system composed of two 
hydrophones, separated in space by some distance, d. 

J 

©«3 

i 2. 

If the output of hydrophone 1 is used as a replica 
with which to compare the output of hydrophone 2, the 
cross-correlation of the two outputs can be obtained. 

What does this do for us? The important property of sea 
noise for this system is that it tends to be isotropic, 
that is to be the same regardless of the direction from 
which it impinges on the hydrophone. In addition, sea 
noise tends to be uncorrelated over distance. On the 
other hand, most signals of interest to us are from a 
single source, usually a submarine, and thus are coherent 
to an extent, and directional in nature. 
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Since sea noise is isotropic, the correlator output 
will be small and relatively constant when no signal 
is present. Now consider what happens in the presence 
o£ a signal o£ interest. The acoustic energy traveling 
outward £rom the source will impinge on the hydrophones 
£rom some given direction. 




1 a 

Following a given wave£ront as it travels across the 
hydrophones, we see that it will encounter one o£ the 
hydrophones be£ore the other unless the direction o£ 
impingement is perpendicular to the line joining the twn 
hydrophones, or the baseline. 

O O 
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The di££erence between the times that the wave£ront 
strikes the two hydrophones depends on the direction 
£rom which it arrives with respect to the baseline. I£ 
the direction o£ arrival is perpendicular to this base- 
line, the wave£ront strikes both hydrophones at the same 
time, and thus there is no time di££erence. I£ the 
direction o£ arrival is parallel to the baseline, the 
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time difference is a maximum. In fact, tlie time difference 
is related to the arrival direction by the following 
equation : 



V cos e ’ 

where v is the speed of sound in the medium, and 6 the 
wavefront arrival angle. 




The cross - corre lation of the two hydrophone outputs 
for a single wave passing the array, shov'^s that the delay 
time T at which a correlation peak appears is a function 
of the wavefront arrival time difference at the hydro- 
phones. If the wavefront arrival angle is 90°, then 
the arrival time difference is zero, and thus the cor- 
relation peak will appear at T = 0. For any other arrival 
angle, there is a positive arrival time difference, At, 
and the correlation peak appears at T = At. This allows 
us to relate the correlation delay time, T; , to the wave 
arrival direction with respect to the array baseline. 

There is a problem of directional ambiguity, as arrivals 
from opposite sides of the array at equal angles with 
respect to the baseline produce the same delay. 
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This ambiguity can be resolved by adding another 
array to cross-fix, or by adding a third hydrophone not 
in line with the first two, to resolve the ambiguity. 




It should be noted that although this method works 
well with broadband radiated signals, it has some draw- 
backs when applied to signals of a single frequency line 
or those with a very broad autocorrelation function. 

The problem is as follows. If a signal consists of a 
single frequency line, its autocorrelation function is 
a sinusoidal function. The question then arises as to 
which peak corresponds to the signal arrival direction. 



141 



aw 

t 



For this reason, a correlation detector is not well suited 
for passive detection of narrow-band signals, but works 
well for broadband signal detection. 

We have now seen how correlation can be of help in 
active systems and in passive broadband signal detection 
systems. What can be done about single- frequency line 
detection? In order to aiiSkvei this question, we examine 
a system which is similar to correlation in many aspects, 
but is not actually correlation. 
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IX. BEAM FORMING 



Beam forming, or electrical steering of arrays, is 
accomplished by a process which is similar to correlation, 
but with one major dif f erence- - it uses summation, vice 
multiplication, of the outputs. We examine the process 
to see how it operates. 

Let us start with a two-hydrophone array. As noted 
in the preceding section, the acoustic wave arrival direction 
determines the wavefront arrival time difference at the 
hydrophones. With this in mind, suppose that it is 
desired to "steer" the array to receive signals from a 
given direction rclativ'c to the array baseline vely . 

that is, to form a "receiving beam" at a given orientation. 



0 

1 



Direction of 
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If we delay the output from hydrophone 2 for a time 

tj = At, and add t, to the time necessary for the acoustic 
d ’ d 

wave from the desired direction to travel from hydrophone 2 
to hydrophone 1, we can selectively bias the array to 
receive signals from that direction. 
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Why does this selectively bias the response of the 
system to signals from that direction? Consider a singl 
impulse arriving from the direction of our selected beam 
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Plotting the output of each hydrophone as a function 
of time produces the following result. 



Hydrophone 2 



Hydrophone 1 




By delaying the output of hydrophone 2 for a time 
t^ = At and adding the output of the two hydrophones, 
the two impulses reinforce each other and the output 
will appear as in the next drawing. 
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If, however, the impulse arrived from a different 
direction, the time difference between arrival at the two 
hydrophones will be different. 
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If we input the output of hydrophone 2 through the 
same delay. At, the plot of the summation of the two 
signals is as shown in the next illustration. 







Thus we see that signals arriving from the selected 
direction (time delay At) will selectively reinforce 
one another in the processor while those from other 
directions will not. 

If we examine a 3-hydrophone linear array we observe 
that ti'ie effect is even liiore pronounced. 




By delaying the output of hydrophone 3 for a length 
of time 2At, and the output of hydrophone 2 for a length 
of time At, a wave from the desired direction will produce 
the effect shown in the next figure. 
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Hydrophone 3 
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It is apparent that adding hydrophones to the array 
makes it more responsive to signals from the chosen 
direction. Also, since the output of the system is not 
a correlation of the individual hydrophone outputs, but 
a simple summation of their energies, the system will 
work for narrow-band or single- frequency line signals 
as well as for broadband signals. 

By using a number of hydrophones in this manner in 
conjunction with an energy detector, and by setting the 
threshold high enough to preclude response to the un- 
reinforced signals from other directions, one has a 
very directional detection capability from a series of 
omnidirectional hydrophones. 
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